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ENGINE MISFIRE DETECTION 

FIELD OF THE INVENTION 
[0001] The present invention relates to internal combustion engines, 
and more particularly to engine misfire detection. 

BACKGROUND OF THE INVENTION 

[0002] Government regulations require automobile manufacturers to 
control the exhaust of engine combustion byproducts such as hydrocarbons, 
carbon monoxide, and nitrous oxide. Emission of such byproducts is typically 
controlled by a catalytic converter. However, if an automobile engine misfires, an 
increased amount of unburned combustion byproducts is passed through the 
catalytic converter. Engine misfire occurs as a result of the absence of spark in a 
cylinder, poor fuel metering, poor compression, or other similar conditions. Over 
time, regular engine misfire can lead to damage of the catalytic converter. 

[0003] A difficulty in the development of an on-board real-time misfire 
detector is the determination of a unified algorithm to detect misfires with high 
detection accuracy and with no false alarms while maintaining a low 
implementation cost, under all required engine operating conditions and misfire 
patterns. These operating conditions include, but are not limited to, varying 
engine speeds, engine loads, road surfaces, gear positions, and mechanical 
dynamics. Furthermore, misfire patterns may be regular, random, or special 
sequences. 
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SUMMARY OF THE INVENTION 

[0004] A method of detecting misfire in an engine comprises detecting 
engine speed fluctuations. A linear model is determined for estimating engine 
firing events based on the engine speed fluctuations. A Kalman filter is applied 
to the linear model to determine parameters of the linear model. A misfire event 
is detected in the engine based on the linear model. 

[0005] In another aspect of the invention, a misfire detection system 
detects misfire in an engine. A sensor determines speed fluctuations of the 
engine. A controller executes a firing event model for estimating engine firing 
events based on the engine speed fluctuations. The controller applies a Kalman 
filter to the model to determine parameters of the model. A misfire detector 
detects a misfire event based on the model. 

[0006] Further areas of applicability of the present invention will 
become apparent from the detailed description provided hereinafter. It should be 
understood that the detailed description and specific examples, while indicating 
the preferred embodiment of the invention, are intended for purposes of 
illustration only and are not intended to limit the scope of the invention. 
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BRIEF DESCRIPTION OF THE DRAWINGS 
[0007] The present invention will become more fully understood from 

the detailed description and the accompanying drawings, wherein: 

[0008] Figure 1 is a functional block diagram of an engine misfire 

detection system according to the present invention; 

[0009] Figure 2 illustrates an engine speed signal in an exemplary six 

cylinder engine; 

[0010] Figure 3 illustrates an engine speed signal with misfire events in 
an exemplary six cylinder engine; and 

[0011] Figure 4 illustrates a signal separation diagram according to the 
present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0012] The following description is merely exemplary in nature and is in 
no way intended to limit the invention, its application, or uses. 

[0013] An engine combustion process is a dynamic system. Normal 
firing and misfire events are inputs that excite the system. Crankshaft speed 
fluctuates in response to the inputs. In other words, crankshaft speed fluctuation 
may be considered an output of the combustion process. Engine misfire 
detection may depend on several vehicle conditions. For example, crankshaft 
speed, average manifold absolute pressure, and firing event signals may be used 
to detect engine misfires. In particular, a current crankshaft speed is related to 
current and previous firing events, as well as previous crankshaft speeds. 
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Therefore, engine misfire detection may be expressed as a general difference 
equation that describes the system in a crankshaft angular sampling domain Ic 
N{k)=H(N,P,l) (equationi) 
Where 

N = [N{k - 1) N{k - 2),...N{k - n)]; 
P = [P{k), P{k - l)...P(Jk - m)]; and 
I_ = [l{k).l{k-\)....l{k-l)] 

[0014] In the above equations, N denotes crankshaft speed, P denotes 
the average manifold absolute pressure (MAP), and I denotes a firing event 
signal. The variables n, m, and / indicate previous samples of the crankshaft 
speed, MAP, and firing event signal, respectively. 

[0015] An inverse model corresponding to the above crankshaft speed 
model may be described in a general fomn as follows: 

l{k)= G(l,N,P) (equation 2) 

Where 

i = [l(k-i),l{k-2),.j{k-l)l 

N = [N{k),N{k-l),...N{k-n)]\ and 

P = [p{k),P{k-\)....P{k-m)] 

[0016] The inverse model implies that the current firing event signal 
may be inversely estimated from the current and past values of the crankshaft 
speed fluctuation, the MAP signal, and previous firing event signals. This is a 
non-linear discrete dynamic system with three inputs and one output. In general. 
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it is very difficult to derive a detailed form of the non-linear dynamic function G. 
Therefore, misfire detection according to the present invention simplifies the 
above inverse model into the following form: 

l{k) = R{n)q{n,p) (equation 3) 

Where 

N = [N{k), N(k - l),...N(k - n)] ; and 
P = [p{k),P{k-l),...P{k-m)] 

[0017] The above model (equation 3) separates the non-linear dynamic 
function G{l,N,P) from equation 2 into r{n) and QiN,P). For the purposes of 
the present invention, the function r{n) is an engine firing event estimator 
function. Misfire signatures are primarily embedded in the function r{n). 
Therefore, the engine firing event estimator function plays a significant role in 
misfire detection. The function Q{n,p) is a load compensator function. 

[0018] The engine firing event estimator function r{n) may be written 
as an [m+if order difference equation as follows: 
y{k) = R{N{k),N{k -l\...N(k - m)) 

= bQN(k)-\-b,N{k - 1)+ ... + b„N{k - m)^v{k) (equation 4) 

[0019] In the above difference equation, bi,i = OX,..,m are the model 
parameters and v{k) represents measurement noise. The y{k) temfi 
corresponds to the l{k) term in equation 3. Furthermore, the Q.{n_,P) function is 
simplified as a unitary gain in the difference equation. Because crankshaft speed 
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iV Is a measured value, the model parameters i),.,i = 0,l,...,w are the only 

unknown variables in the difference equation. Various system identification 
techniques may be used to estimate the model parameters of the difference 
equation. The engine misfire detection system of the present invention uses a 
Kalman filter to estimate the model parameters. 

[0020] To use a Kalman filter to estimate the model parameters, the 
difference equation (equation 4) is first reformulated using standard state space 
systems equations as shown below: 

x{k + 1) = F{k)x{k) + B{k)uik) + G{k)wik) (equation 5) 

yik) = H(k)x{k) + v(k) (equation 6) 

Where 

x^ik) = bo{k) + w^ik) 
X2(k) = bi(k) + W2ik) 

Xn„iik) = bjk) + w„,,(k) 

[0021] In the above equations, w^ik) is a zero mean Gaussian noise 
vector, wherein /= 1,2,3,. ..,mf1. Therefore, the result is jc(ii:+l) = ;c(jk) + w(jt), 
where 

x{k) = (equation 7) 
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[0022] The result of the above equation is a linear system equation that 
is equivalent to equation 5 with Fik) = G{k) = / , where / is a unitary matrix and 

M(ik) = 0. 

[0023] Referring again to equation 6, if 

Hik) = [Nik) N(k-l) Nik-2)...Nik-m)], then equation 6 represents the 
same model as equation 4. Ultimately, the present model for misfire detection is 
represented by equations 6 and 7. In other words, equations 6 and 7 represent a 
firing events model. A Kalman filter estimator can be derived based on the 
model. The base equations of a Kalman filter are: 
xik + l/k) 

= Fik)x{k /k-l) + B(k)uik) + G(k)wik) 

+ [Fik)I,ik/k- 1)H (ik) + G(/k)5 (k)] [H (*) Z()k / ik - \)H ^ik) + Rik)]'' 
[y(k) - H(k)x(k /k-l)- v (/k)] 
Jc(0/-l) = 3c(0) 

[0024] Where K{k) is the Kalman gain given by: 
Kik) = [Fik) Z(ik / ik - (ik) + Gik)Sik)][H ik) / ik - l)H^ (k) + Rik)]-' 

[0025] The base covariance equations are: 
I(ik + l/ik) 

= Fik)I.{k/k-l)F^{k) + Gik)Qik)G^ {k) - K{k)[Fik)I,(k/ k-l)H^ {k) + G{k)S{k)Y 
E(0/-1) = Z(0) 

[0026] The noise statistics assumptions are: 
E[wik)] = wik) 
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E[v{k)] = m 

E{ [w{k) - - vv(;)f } = QOOd^ > 0 

E{ [vik) - v{k)mj) - vu)f } = ms^ > 0 

E{ [wik) - w(*)][v(;) - v( } = S{k)S» ^ 0 
£{[;((0)-3c(0)][v(/:)-v(^)f } = 0 

[0027] Using the systems equations 6 and 7, a predictive Kalman filter 
equation can be simplified to: 

x{k + \lk) = x{klk-\) + K{k){y{k) - H{k)x{k Ik -I)] (equation 8) 

a;(0/-1) = x(0) 

I(/t + 1/ Jt) = [/ - K{k)Hik)] Y.{klk-\) + Q{k) 
Z(0/-1) = S(0) 

[0028] The vector x{k) in equations 6 and 7 can be estimated from 
actual vehicle test data. This vector determines the "tf coefficients of the firing 
events model. The H{k) term can be defined based on engine speed N. 
Therefore, the firing event signal yik) from equation 6 can be determined based 
on crankshaft speed N and the application of a Kalman filter to estimate the 
vector x{k) . 

[0029] A misfire detection system 10 according to the present invention 
is shown in Figure 1. An RPM sensor 12 senses a speed of an engine 14. The 
RPM sensor generates an engine speed signal and inputs the engine speed 
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Signal to a signal processing system 16. The signal processing system 16 
processes the engine speed signal and derives a firing event signal y(k) 18 
according to the above equations and known engine characteristics. For 
example, the signal processing system correlates the engine speed signal to 
crankshaft speed. The firing event signal 18 is input to a Kalman filter 20. 

[0030] The Kalman filter 20 estimates unknown parameters of the firing 
event signal 18. A misfire detector 22 determines engine misfire events 
according to an output of the Kalman filter 20. Additionally, the misfire detector 
22 may consider other engine conditions to compensate for signal output 
variations caused by factors such as engine load. For example, an average 
manifold absolute pressure and the engine speed determine the Qiji^P) term as 
shown in equation 3. The misfire detection system 10 may use a lookup table or 
a surface mapping function to determine Q{fi,P). Alternatively, other suitable 
methods may be used to determine gQv,p). For example, the signal processing 

system 16 and the misfire detector 22 may be combined in a single device. 

[0031] Figure 2 is an example of a six cylinder engine at 2000 RPM 
with no load and no misfire events. An engine speed signal depicts a plot of 
RPMs against a firing number. Figure 3 illustrates a six cylinder engine speed 
signal with no load and 20% cylinder misfire. Figure 4 is a signal separation 
diagram of the engine speed signal resulting from the present invention. The 
signal separation diagram plots misfire events and normal firing events against 
firing number, where "O" denotes a normal firing event and "X" denotes a misfire 
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event. This plot depicts a signal to noise ratio of the misfire detection system 
according to the present invention. 

[0032] The description of the invention is merely exemplary in nature 
and, thus, variations that do not depart from the gist of the invention are intended 
to be within the scope of the invention. Such variations are not to be regarded as 
a departure from the spirit and scope of the invention. 
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